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Abstract 

The Ul{ 3) x Ur{ 3) Linear Sigma Model (LSM) with quark degrees of freedom is 
used to show that radiative corrections generate undetermined finite contributions. 
Their origin is related to surface terms which are differences between divergent 
integrals with the same degree of divergence. The technique used to detect these 
ambiguities is an implicit regularization on basic divergent integrals that do not 
depend on external momenta. We show that such contributions are absorbed by 
renormalization or fixed by symmetry requirements. The general expression for 
surface terms is derived. Renormalization group coefficients are calculated, as well 
as relevant observables for this model, such as ff. and the pion and kaon form 
factors. 
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I. Introduction 


In dealing with ultraviolet infinities which stem from the semi-classical 
expansion in a quantum field theory, the renormalization program plays an 
essential role to fit the experimental data which the model is meant to de¬ 
scribe. In regularizing and subtracting a divergence there appears a hidden 
parameter which plays the role of renormalization group scale in renormaliz- 
able models. The imposition that the Green functions do not depend on how 
the infinities are subtracted from the amplitudes can be formulated by de¬ 
manding that the Green functions are independent of such parameter. This 
gives rise to the running couplings, quantum mechanical mass generation, 
etc. 

Recently Jackiw [T] raised an interesting matter concerning arbitrary pa¬ 
rameters which appear in perturbative calculations in certain field theoret¬ 
ical models. Such arbitrary (regularization dependent) parameters show up 
as differences between divergent integrals of the same superficial degree of 
divergence. Usually physicists do not like to have free parameters in their 
theoretical models. The prescription to deal with such arbitrary quantities 
is threefold: 1) check whether they may be considered as finite counterterms 
(as it is usual in renormalizable models); 2) check whether symmetry con¬ 
straints (e.g. supersymmetry, gauge symmetry, etc.) may require that such 
parameters assume a definite value and 3) see if as a genuine free parameter, 
phenomenology can make use of it to fit the data in the underlying model. 
Of course such parameters are expected to play an essential role in finite 
theories and effective models. 

A well known example where a finite ambiguity plays an essential role 
is in the description of the Adler-Bardeen-Bcll-Jackiw anomaly |2], |3j- In 
this case the finite constant is a manifestation of an important symmetry 
breaking at the level of perturbative quantum corrections. A democratic 
description of the anomaly between the vector and axial sectors is achieved 
when we allow a regularization dependent parameter to be arbitrary [1]. 

The ideal arena to implement such ideas in the momentum space is the 
Implicit Regularization (IR) scheme [21 , which has already been successfully 
applied to several examples 0 , 0 . 0 . 0 , 13 ]. The idea behind IR is that one 

(a) separates external momentum independent divergencies from finite 
contributions, without the need of an explicit regulator; 

(b) leaves their divergent content in the form of basic integrals until the 
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end of the calculation; 

(c) does not modify the underlying theory, preserving its symmetries and 
space-time dimension. 

The purpose of the present contribution is the investigation in the LSM 
PI of the role played by ambiguities in the context of the above discussed 
ideas. We show that radiative corrections yield together with well known 
infinities also finite undetermined contributions. In this work we will follow 
Jackiw’s suggestion and leave the arbitrary quantities open until the end of 
the calculations. Then, we will show that these undetermined contributions 
can be fixed by symmetry relations or absorbed by renormalization, in agree¬ 
ment with the ideas presented in m- From the point of view of the current 
phenomenology, because the present Lagrangian does not contain the impor¬ 
tant ’t Hooft determinantal interaction {T2], we restrict the fits to the pion 
and kaon sectors, which are not crucially dependent on it. 

The work is organized as follows: in section 2. the method is presented for 
the case of amplitudes containing different masses in the fermion propagators 
and we show the origin of ambiguities. The renormalization group scale is 
also introduced within IR. In section 3. the model is renormalized and some 
ambiguities are shown to be absorbed by renormalization. Moreover the (3- 
function is explicitly calculated and the renormalization group equation for 
the pion coupling constant is solved in section 3A. The position of the well- 
known Landau pole is identified and coincides with the position of the pole 
in the large euclidean momentum region of the pion propagator. In section 
4. we derive the expressions for the weak decay constants f n and / K , as well 
as the electromagnetic form factor for the kaon. Numerical results are given 
in section 5., and conclusions are given in section 6. 


2. Treatment of divergent integrals, ambiguities and renormal¬ 
ization group scale 


In the present section we illustrate the relevant technical details of the 
renormalization procedure we use by working out explicitly a one-loop Feyn¬ 
man amplitude. 
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Consider the pseudoscalar amplitude with two quark masses m u and m s 

nP? = 1 L |y Tr ' (1) 

where the symbol A stands for a regulator which needs not be explicitated, 
but is necessary to give a meaning to n pp . One is then allowed to alge¬ 
braically manipulate the integrand. We do it in such a way that the diver¬ 
gencies appear as integrals and are separated from the finite (external mo¬ 
mentum dependent) contribution to Eq. After taking the Dirac trace, 
we use the following algebraic identity 


1 = 1 _ P 2 ~ 2p-k _ / 

(k — p) 2 — m 2 k 2 — m? [(p — k) 2 — m 2 ](k 2 — m 2 ) 

at the level of the integrand. Note that it allows one to confine the external 
momentum dependence in more convergent integrals. This relation should be 
used recursively until the finite part is completely separated from divergent 
integrals. We get 

PP = f f d A k i r d A k i 

\J (2vr Yk 2 -m 2 u J (2vr ) A {k - p) 2 - m 2 

+ [K - m u f - p 2 } J ( fc 2_ m 2)[( fc _ p )2_ m 2] } • ( 3 ) 


The first integral on the RHS is what we call a basic quadratic divergence 



d A k i 
(27r) 4 k 2 — m 2 


( 4 ) 


The second integral on the RHS is also a quadratic divergence, but it still 
possesses an external momentum dependence. If one uses Eq. (( 21 ), one sees 
that an arbitrariness emerges 


T(p 2 , nrii) = j 


d A k i 
(27r) 4 (k — p) 2 — m 2 


I q{m 2 s )+p^p v A^, 


where 




d A k Akf,k u . r d A k g 
( 2 tt ) 4 (. k 2 - mlf 1 J ( 2 tt ) 4 (, k 2 - m 2 ) 2 ~ 


( 5 ) 

( 6 ) 
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and a = const, i.e. the difference between two logarithmically divergent 
integrals. Eq. © is a surface term, a particular case of the general expression 
for surface terms (see derivation in Appendix ) 


fj f A t lA t 2---M2n _ ^ r(A 1) P[liPl_l2 ■ • ‘P^2n 

i 2 — m 2 ) A ~ n ~ l T(A -n-1) (p 2 — m 2 )" 4-1 

(7) 

where u is a continuous dimension, /ij is a Lorentz index and <7 Ml/i2 ... M2n is the 
standard notation for symmetrized products of the g, w tensor. It is clear that 
at u = 2, A = 4 and n = 1 we reproduce Eq. ©. In dimensional regular¬ 
ization or Pauli-Villars, following the usual prescriptions to preserve gauge 
invariance, one would get zero for However for other regularization 

prescriptions one may expect a dimensionless number. Following Jackiw’s 
suggestion we leave this arbitrary number open until the end of the calcula¬ 
tions and only then see if it must be fixed by Ward identities. In QED and 
QCD gauge invariance forces it to be zero Hj- In the anomalous pion decay, 
on the other hand, it must be finite [1]. 

The third term in Eq. © is also logarithmically divergent. In the spirit 
of Implicit Regularization the integral can be displayed as 


Cpt-2) = . /' 
AtlM2..-M2n / 


d 2u p 

Wf 


d 4 k 


(27t) 4 ( k 2 — m 2 ) \{k — p) 2 — m 2 


d A k 


(27r) 4 ( k 2 — m 2 ) 2 
■ d 4 k i 

(27t) 4 ( k 2 — m 2 ) 2 
• d 4 k i 

(2tt) 4 (k 2 - C 2 ) 2 


i rV 2 2 2 2\ 

■ + F(m u ,m s ,p ;mj 

+ F(m 2 u ,m 2 s ,p 2 -,m 2 s ) 

+ F(m 2 u ,m 2 ,p 2 - C 2 ) 
( 8 ) 


where again we separated another basic divergent integral 

d A k i 


hogim 2 ) = j 

and a finite contribution 


(2t r) 4 (k 2 - m 2 ) 2 ’ 


(9) 


C 2 ) = ^ [' (10 ) 


(47t) 2 Jo 


< 2 


In this expression the last argument corresponds to a squared mass in the 
denominator of the logarithm and plays the role of a scale. This comes from 
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the property of the logarithmic divergence, 


1 771 ^ 

Wm 2 ) = «C 2 ) + ^M^), (11) 

and its intimate connection with the finite part of the amplitude. The arbi¬ 
trary scale Q 2 represents in this method the renormalization group scale. 
Finally we obtain that 

n PP = 2 {I q (ml) + I q (m 2 s )+p 2 a 

+ [{m u -m s ) 2 -p 2 ](I log {C) + (12) 



(13) 


(14) 


3. The model 

We start by considering the following generating functional 

Z = I ]\Vcr 0a Vir 0a VqVq exp(iS(q,q,a 0 ,n 0 )) (15) 

** a 

with the action S(q, q, cr 0 , 7To) of a £7(3) x U( 3) linear sigma-model including 
fermionic degrees of freedom 

S(q,q,<Jo,n 0 ) = J d 4 x C q - ^ tr [(BB ] ) 2 } + tr(c 0 a 0 ) 

+ y tr(d^cr 0 d^a 0 + d^od^no) , (16) 

C q = q[i ^ - g 0 (a 0 + i^f 5 7r 0 )\q. (17) 
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Here B = a 0 + *7r 0 , B~ f — cr 0 — Z7r 0 , used in the definition of the quadratic and 
quartic invariants in mesonic fields m 

tr {BB ] ) = tr(crQ + 7 Tq) 

tr l(BB ] f] = tr{(crQ + tt ^) 2 - [cx 0 , tt 0 ] 2 }, (18) 

with the scalar and pseudoscalar fields cr 0 = A a cxQ and 7 r 0 = A a 7 TQ containing 
the standard 17(3) matrices A a (a=0,1...8) in flavor space. The subscript 0 
stands for bare quantities, C q describes the coupling of quark fields g, q to 
the mesons, go is dimensionless. 

The symmetry transformation properties of the fermionic fields are 

5q = i( a ' + 7 5 (3')q, 

8q = —iq(a' - y 5 f3') (19) 

and for the mesonic fields 

5ao = i[a', cr 0 ] + {/?', 7T 0 } 

Sttq = i[a ', 7T 0 ] - {/?', cr 0 } (20) 

where the parameters of the infinitesimal global transformations a' and (3' 
are Hermitian flavour matrices. The explicit symmetry breaking piece is 
introduced via the linear term proportional to Cocro with a diagonal matrix¬ 
valued strength Co = diag{co u ,co u ,cos}- Without the quartic interaction 
and kinetic terms for the mesons, proportional to /q and Ao g respectively, 
the action corresponds to a Nambu-Jona-Lasinio (NJL) type lagrangian [TT] 
in semi-bosonized form; this is achieved by the method of auxiliary fields 
na. which has been widely used in various extensions of the original NJL 
lagrangian, see e.g. [2DH22|. We relate the cq to the diagonal current quark 
mass matrix m 0 through the choice Co = g^rho- We restrict our study to the 
case rh u = rh^ ^ m s , which breaks the unitary SU(3) symmetry down to the 
subgroup 517(2)7 x U{1)y (isospin-hypercharge). Upon the exact integration 
over the quark variables in the path integral Eq. m one obtains for the real 
part of the effective action associated with the fermions (we do not consider 
here anomalous processes), 

J VqVq exp(i J d 4 xC q ) —■> exp(ln | det D E \) = exp ^-Tr Ih^D^De)^ ( 21 ) 
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where we use the strictly positive unbounded Hermitian second order elliptic 
operator D^ E D E , which is chiral and gauge covariant na-na and where Tr 
designates functional trace, D E stands for the euclidean Dirac operator 

D e = - g 0 (a 0 + Wo), 

d e = -*7#A “ 9o(cr 0 - i 7 5 vr 0 ). ( 22 ) 

Since the scalar fields possess non-vanishing vacuum expectation values, we 
perform the shift <7 0 —> ^ + ctq, where m represents the hnite constituent 
quark mass matrix and obtain after this shift 

D ] e D e = m 2 - d 2 + Y, (23) 

with the background mesonic helds Pi 


777 / 777 / 

Y = + *75^0) + 9 o ( o-q + (—1 °o} + ^0 + *75 [00 H-, ttq] ) , 


9 0 


9 0 


(24) 


leading to the expansion 


Tr(ln D%D e ) = -Tr{(-<9 2 + m 2 )[l + (- d 2 + m 2 ) 1 [^o 7 M ( 9 A 1 ( 7 o + * 75 ^vr 0 ) 


n=l 


n 

2 1 2 


+fi , o( CT o + {cr 0 , m/go} + 7r 0 + *75^0, 0o + m/ 00 ])]]}" (25) 


From now on it is convenient to use the following representation of the 
scalar and pseudoscalar helds 



/ 

(7 u 

V2 

a + 

< 5 + \ 


a 


5 ° 

V 

5~ 

5 ° 

£s_ 

V2 ) 


~7T 


f 

v7 

7T + 

iF+ 

\ 


71 

4>d 

W° 


V 

K~ 

V " 

J?° 

< i> s 
V2 

/ 


(26) 


The gap equation will be obtained by considering the n — 1 term in the 
expansion Eq. (EH) and the linear contributions in 07 from the remaining 
terms of Eq. (d after the mass shift. We get 


ho 

9 o 


Coi 


m * ~ "T - 8 N c mil q (mi) + 2 = 0 , 


9 o 


9 0 


where {i = u,s}. 


(27) 








The n — 2 term in the expansion contains all the other divergent con¬ 
tributions which go up to four-point functions. In the calculations there 
appear the two basic divergent integrals I q and Ii og , Eqs. (jH) and 0, as 
well as the difference between two logarithmically divergent integrals, a, Eq. 
Q . Expressions m and after properly taking into account go and 

trace factors in color (N c ) and Dirac spaces, are the amplitudes, continued 
to Minkowski space, resulting from terms quadratic in the fields of 


In the following we need the wave function renormalizations, which are 
readily obtained as the coefficients of the p 2 terms in all expressions quadratic 
in the fields, stemming from the kinetic terms proportional to ff and from 
the amplitudes (FT21) . (fni) 

Z~ x (m h C 2 ) = fo ~ 4 N c gl{Ii og (C ) ~ a + F(m;, rrij, 0; C 2 ) 

- (rrii +mj^F^m^rrij, 0;C 2 )}, (28) 




fo ~ 4N c gl{I log (C) - a + F(m u m j: 0; C 2 ) 

(mi - mjfF'frrii , rrij,0 ; C 2 )}- (29) 


where F' represents the derivative with respect to square external momen¬ 
tum of F, taken at p 2 = 0. The retained finite terms F(rrii,mj, 0; (y) and 
F'(mi,mj, 0; C 2 ) lead to different wave function renormalizations for the dif¬ 
ferent members of the pseudoscalar and scalar nonets; the higher contribu¬ 
tions to the expansion of F(rrii, rnj,p 2 ; £ 2 ) will be absorbed in the renormal¬ 
ized propagators of the mesons (see Eq. ipp below). 

The renormalized coupling constants g a (rrii,mj) and rrij) are then 

obtained as 


9o 


glfm^rrij-X 2 ) 


= Z a (m^mfX ) 


(30) 


9l 


9n( m i, m j'i C 2 ) 


= z n (rrii, rrij] ( ), 


(31) 


and the renormalized masses g a ^(rni,mj]( 2 ) become 


Z a , v (mi, rrij ; C 2 ) 


Ma,w(mi, my, C 2 ) 

t4 ~ 4N c gl((I q (m 2 ) + I q (m))) 
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- 4 N c g%(mi ± m j ) 2 (/ /og (( 2 ) + F(mi, m,, 0; C 2 )) 

+ 2^(m 2 + m 2 ± m, ; mj), (32) 

9o 

where M 2 goes with the plus signs, and the renormalized quartic coupling 


X q _ Xoq 

gl(m u ,m u ]C 2 ) g$ 


-m c i log {c 2 ). 


(33) 


We have retained in the renormalization of the masses the finite terms 
F(mi, rrij, 0; ( 2 ). Note that all ambiguities a appear in the held renormal¬ 
ization coefficients (I28H29D . In principle they could assume different values 
for the different processes. Chiral symmetry restricts them to have the same 
value, as we shall see below. Using relations CD and (El one obtains that 
the kaon and pion renormalization constants are related as (where we have 
identified g n = g^{m u , rn u : ( 2 ) and g K = g n (m u , m s \ ( 2 )) 


— = — - 4 N c [F(m u , m s , 0; ( 2 ) - F(m u , m u , 0; ( 2 ) 

gi gi 

- (m u -m s ) 2 F'(m u ,m s ,0-,( 2 )\ 

Xn 


(4 h 2 

9k 


= -f + 2 m s (m s - m u ) 
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9, 


- 4:N c F(m u ,m s ,0-X ) 


F(m u , m s , 0; £ 2 ) = 


(4tt) 


m~ 


rn'i — mi 


In 


m~ 


m: 


+ In 


mz 


With these definitions we obtain the effective action 


(34) 

(35) 

(36) 


S 


Skin “1“ S mass “I - 


inti 


(37) 


where S^n , S mass , S int are the actions related to kinetic, mass and interaction 
terms of the mesons (in the equations for Skin, S mass below we omit the index 
0 in the bare fields in order not to clutter the notation). 

Skin = I d A x rn u ) {d IL a u d ,l a u + d^a d d^a d + d ll <j + d l ‘a~) 

+ Z~ 1 {rn s ,m s )d lx a s d tx a s 

+ Z-\m u ,m s ){d^ + d^~ + c^W) 

+ Z~ 1 {m u , m u )(a M 0 u a /i 0 u + dpfadPfa + ^vr+^Tr - ) 

+ Z~ 1 {m s ,m s )d^ s d tl ^ s 

+ Z;\m u , m s )(d lx K + d tl K~ + d IJ K () d"K Q )} , (38) 
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Smass = - J d A x {M 2 (m u , m u )(crl + cr} + a + a ) + M 2 (m s , m s )a 2 
+ M 2 (m u , m s )(5° 5°+ 5+5~) 

+ M 2 (m u , m u )((j) 2 + 0^ + 7r + 7T“) + Ml(m S} m s )(t) 2 s 
+ M 2 (m u , m s )(K°K° + K+K~ )} . (39) 

If we use the gap equations in order to eliminate the quadratic divergencies 
in the expressions for the renormalized masses /^^(m*, m^), Eq. (El , it is 
possible to define the renormalized coupling C t through 


C u 

C s 


9l 

CsO 

9o 


= m „ 


= m,. 




2 /£ h 2 


9l 


9 , 


2 (m t 


"u) 2 (^§ - 4-/V c F(m u , m s , 0; ( 2 )) 
£ 2 


(40) 


Using El one obtains also 


C u + Cs — (m u + m s )^-F. (41) 

9 k 

Now we write out the renormalized propagators for pions and kaons, since 
they will be used in the remaining, 


A n \p 2 ) = p 2 ~i4 + 4:N c glF fin (m u ,m u ,p 2 -,( 2 ) (42) 

A^(p 2 ) = P 2 ~ pl + 4:N c g 2 j: fin (m u ,m s ,p 2 ; C 2 ) (43) 


with the hnite momentum dependent contributions 

Ffi n (m u , 7n u ,p 2 ] ( 2 ) = -p 2 (F(m u ,m u ,p 2 ]( 2 ) - F(m u , m u , 0; C 2 )) (44) 

where in particular F(rrii, rrii, 0; = 0, 

^f in (m u ,m s ,p 2 ; ( 2 ) = {(m s -m u ) 2 - p 2 }{F(m u ,m s ,p 2 ]( 2 ) - F(m u , m s , 0; C 2 )} 
- p 2 (m s - m u ) 2 F'(m u , m s , 0; ( 2 ) (45) 
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where use has been made of the normalization conditions 


A„( 0 ) = -It 


2 

7 


dA~i 


0P 2 


dp 2 


=n — 1- 


Ip —0 


( 46 ) 


We obtain finally the physical pseudoscalar masses as zeros of these propa¬ 
gators 


m l = l4- 4N c9lFfin(m u ,m u ,ml;C 2 ) 

m l = 9l~4N c g 2 K J: fin (m u ,m s ,ml;( 2 ) (47) 


At the physical meson masses we obtain the following pseudoscalar quark 
couplings 


dF fin (m u ,m s ,p 2 ]C 2 ^ 
dp 2 l? 

dYt fin (m u ,m s ,p 2 ] C 2 )] 
dp 2 


-2 


9 nan ~ 9 -k +^N C 


iirqq 

— 2 


9 nqq ~ 9 k + ^ N c 

3.A - The (3 -function 


(48) 


In order to illustrate the role of the arbitrary scale ( 2 introduced in the 
previous sections, we evaluate the /3-function for the pseudoscalars 

A = cJtStK 2 )- (49) 

Using the Eq. CO we get to one loop order 


Ar = 



(50) 


,/From the renormalization group equation PI) and © one solves for g n 
and gets 


9i(n = 


gi( e 


£^(C 2 )/n 


- £/2 
A 


(51) 


Analogously to QED at one loop level, a Landau pole appears in equation 
© and it corresponds to a pole in the large Euclidean momentum region 
of the pion propagator. This is easily seen by expanding for large negative 
p 2 the inverse of the renormalized pion propagator, Eq. (S3) 


A 


-i 


i/k * 


4(/r 2 




(52) 
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4. Coupling to external fields 

In order to make contact to phenomenology an important issue is the 
calculation of electromagnetic form factors and weak decay constants. We 
treat the vector and axial currents as external classical fields. We extract 
these quantities from a generalized expression for the Dirac operator, now 
containing the vector and axial-vector A^ fields. The previous expression 
for D^D is now replaced by 


D ] e D e = m 2 -dl + igo-f^A^ao - i'y sA^ttq) - ^[7^, qjT^ 

+ fi'o((°o + Wo,m/go} + 'Ko+i 7 5 [a 0 + m/g 0 ,'K 0 }), (53) 


where 

- dj = -a; + 2 tr„a„ + + rj ( 54 ) 

and T m = + 75 A^, F^ = + 75 A^ v . The covariant derivatives are 

V M cr 0 = d^a 0 + {A e , vt 0 } - i\Vn, a 0 + m/g Q \ (55) 

V m 7t 0 = <9^77) - {A^, a 0 + m/go} - i[V^, ^ r 0 ] (56) 

= d^Vu ~ d v V» - i[A„, A v \ - i[ V„, K] (57) 

A,u, = d E A„ - - i[V iu , A v \ - i[A E , V v ). (58) 


The calculation of f n and f K to one loop order is now straightforward, 
after collecting the relevant vertices involving one pseudoscalar and one axial- 
vector fields. Since we are treating the external currents classically, we disre¬ 
gard all quantum dynamical corrections to the vector fields. Also the kinetic 
terms proportional to /g involving the scalar and pseudoscalar fields in Eq. 
m have to be minimally coupled to the external currents through the co¬ 
variant derivatives Eq. We obtain the bare /° as 




/;K) = —(-AN c g 2 I log (C 2 ) + f 2 -AN c g 2 F(m u ,7n U: mi,C 2 )) 

9 0 

Ffin(m u ,m u ,ml,( 2 ) 


m. 


U Z n 1 (m u , m u ) I 1 - 4 N c g' 2 (a - 


go 


rnt 


(59) 
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In order to obtain the Goldberger-Treiman relation of chiral symmetry 
it is necessary to avoid the ambiguous a term present in Eq. ©. This can 
be done by either choosing a = 0 at this point j‘23j or by introducing the 
following counter-terms 

f 2^2 rryi 

Lc.t. = -^■Tr(-d ll Tr{A IJl ,a 0 -\ - } - ^vr 0 [V^, tt ] (60) 

4 9o 

TTl 

+dfj,(Jo{An, 7To} - id^a 0 [V^ cr 0 H-]), 

go 

with a parameter / 2 to be fixed appropriately. This is the minimal chiral 
and gauge invariant combination available which contains a term that can 
absorb the a ambiguity. For the present purpose this is enough, although 
one might need to consider higher than cubic field interaction terms as well 
to realize the complete renormalization program. Note that this counterterm 
has consequences for the scalar-vector (SV) sector, which will be discussed 
below. Including the terms from © that contribute to the pion weak decay 
constant (proportional to A^d^n), we obtain 


fw( m n) = —Z-\m u ,m u ) (l + fgl-*N c gl(a 

9 o V 

Using that 


F fin (m u ,m u ,ml,C z 

mi 



4 N c gl 


mi 


fin\ 


:,c 2 ) = i- 




mi 


(62) 


and the relation between the physical and bare pion field at the physical pion 
mass, 7r = 7 r p,, 2 zaa ; we g e ^ ^} ie renormalized weak pion decay constant 


2 

r t 2 \ g-wqq AG 

MmJ = m u —^— 


qt mi 

«^ 7 T 71 


(63) 


provided 

f = 4 N c a (64) 

in Eq. m, as a consequence of requiring chiral symmetry. 

Now, for the kaon weak decay constant one obtains 
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- ^N c g^F(m u ,7n s ,m 2 K ;C 2 ) , 


The term proportional to (rn u + m s ) is the term which in the limit of equal 
masses reproduces f n . We see that the introduction of the counter-terms (EOD 
in the model Lagrangian, with p = 4 N c a allows the (rn u + m s ) part of (E3) 
to be renormalized in the same way we did for f n . However there remains 
an a dependent piece, proportional to the difference of the quark masses, 
which cannot be removed by the counterterm, nor by renormalization. It is 
interesting to note that the ambiguity a can neither be fixed by the Ward 
identity related to the pseudoscalar/axialvector radiative amplitude 



where 



(67) 


and T(p 2 , mi) given in Eq. One can verify that Eci. flbbl) is fulhlled with¬ 
out fixing ct, as long as the ambiguity present in the pseudoscalar/axialvector 
amplitude and the one present in the pseudoscalar/pseudoscalar amplitude 
are the same. One concludes that in the PA sector alone an ambiguity 
remains undetermined, which is compatible with the Goldberger-Treiman 
relation and modifies the PCAC relation only as a next to leading order cor¬ 
rection in the symmetry breaking term and is leading correction in flavor 
symmetry breaking. This conclusion is however fallacious, as the chiral sym¬ 
metry partner sector SV has to be analyzed in parallel. By fixing a through 
Eq. m, Eq.® introduces counterterms in the scalar-vector sector 



where V+ = V 4/x - iV 5 ^, 1/ = + iV 5/l , V'p = V 6fl - iV 7/l and = 

W/j, + fVV/j,- These counterterms absorb the ambiguous a proportional terms 
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(69) 


which are generated radiatively through 


n 


fi,ab _ 
SV ~ 


i 

4 


d 4 k 


Tr 


-7^A 


fc — m /t — jj) 


m 



However, opposite to the PA case, these counterterms are proportional 
to the difference of the quark masses and will leave untouched an ambiguous 
term proportional to the sum of quark masses. This latter appears in a 
similar fashion as in Eq. (ESI) through the difference of quadratic divergencies 
with different external momentum dependence. This term violates chiral 
symmetry, since the SV amplitude must vanish in the limit of equal quark 
masses. We conclude therefore that the ambiguity a must be set to zero 
in this sector. This immediately implies that it must be zero as well in the 
PA sector, as the counterterms in the two sectors are linked through chiral 
symmetry. 

We consider from now on the case with a = 0. 

Returning to Eq. ESI), we observe that the quadratic divergencies can be 
removed through the gap equations and the other divergencies by renormal¬ 
ization. Using Eqs. EHESI) and (ED f K reduces to 


hirnl) 


9nqq~ 


{m u + m s ) 




m l9l 


(70) 


At this stage it is useful to see how the expressions for and f K comply 
with PCAC. For that we use the explicit symmetry breaking term of the 
Lagrangian in Eq. (USD, m 


&£sb = —STr(c 0 a 0 ) 

2go 

to obtain for instance for one isotopic component of the pion held 

~TuT~ = -27ri TT = g ^ Cu 

up i 9o 

and for one flavour component of the kaon 

dSCsB Cou + Cqs _ ^ph 


—— = -Kt- 

ofl 4 90 


= gKqq{C u + C.) 


(71) 


(72) 


(73) 


Upon using here Eqs. (ED and ED and the expressions El and (ED one 
obtains 

dS Q^ B = -2 7rf ml f n (ml) (74) 
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Set 

m u = 2ltiMeV 

m u = 350 MeV 

9n 

9 71 

m s 

fn 

< rl > 1 / 2 

2.250 

Ul.lMeV 

479.0 MeV 
lU.OMeV 
0.573 fm 

3.752 

Ul.OMeV 

601.7 MeV 
71.6 MeV 
0.608 fm 


Tabic 1: Parameters of the model, fitted to reproduce m n = 139 MeV, m K = 
494 MeV and f n = 93.3 MeV , with m u = 210 MeV and m u = 350 MeV. 

= -2 Kfm 2 J K (ml). (75) 

So, the pion and the kaon fulfill exactly PCAC. 


5. Numerical Results 

In this section, we fit the model parameters to reproduce the pion and 
kaon masses, as well as the pion weak decay constant, in order to obtain 
a numerical estimative to the kaon weak decay constant and the kaon elec¬ 
tromagnetic form factor. The model parameters are fixed by means of Eq. 
m and (0) in order to reproduce m n = 139 MeV, m K = 494:MeV and 
f n = 93.3 MeV, with the renormalization constants g\ and fi 2 given by (ED 
and (1351) . 

As pointed in ref. [23j, the action defined by (EZD gives the same formal 
results as the conventional NJL model (with a A —» oo cut-off) in the case 
where X q — 1. In particular, the Nambn relation in the chiral limit rn a = 2 m 
holds only in the \ q = 1 case. We will proceed the parameters fitting with 
X q = 1, and will return to this point latter. We also use m u = 210 MeV and 
m u = 350 MeV , values largely employed on literature, generating the two 
set of parameters showed on table I. 

With these choices for the parameters, and with N c = 3, we obtain for 
the kaon weak decay constant, f K = 114.2MeE (for m u = 210MeE) and 
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f K = 71.6 MeV (for rn u = 350 MeV). 

The value of f K estimated by the set of parameters corresponding to 
rn u = 210MeE is only ~ 3% lower than the experimental results. Neverthe¬ 
less, for m u = 350 MeV the kaon weak decay constant is still away from its 
experimental value, as already known. 

We also evaluate the pion and kaon electromagnetic form factors, and the 
corresponding pion and kaon charge radius in the space-like region, given by 


F k = ~N c g Kqq [— + F(ml , m 2 u , q 2 ; C ) - (m 2 K - (m s - m u f)F'(m s , m u , m- Q 2 ) + 
9tc 


It(q , m s , m u )] - N c g Kqq [— + F(m s ,m s , q ; C ) - 

9n 

(ml - (m s - m u ) 2 )F'(m u , m S) m 2 K \ C 2 ) + It(q 2 ,m u ,m s )}, 


(76) 


and 


where 

I t (q 2 ,m u ,m s ) 

with 


< r; >= —6 


dF K (Q 2 ) 
dQ 


— \Q 2 ^0 


d A k 


(27r) 4 ( k 2 — rn 2 )[(k — ki) 2 — m 2 \[(k — k 2 ) 2 — m 


(77) 

(78) 


k\ = k 2 2 = ml (79) 

and 

k l .k 2 = ml~ < ^. (80) 

Also, in Eq. (EH) and in figures 1-3, Q 2 = —q 2 as usual. 

For the pion electromagnetic form factor and pion radius, the expressions 
are given by dZHD and m in the m s = m u limit. 

In figure 1 we plot the kaon electromagnetic form factor in the space-like 
region, evaluated in the present model compared with experimental data. 
This result shows a good agreement between the model and experiment. In 
figure 2 we focused the region near Q 2 = 0, that highlights the differences in 
the slopes of the m u = 210 MeV and m u = 350 MeV fittings, and explains 
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Q 2 (MeV 2 ) 


Figure 1: Kaon electromagnetic form factor in the space-like region com¬ 
pared with experimental data j22J • Solid line corresponds to the set of param¬ 
eters fitted with rn u = 350 MeV, dashed line corresponds to m u = 210MeV. 
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Q 2 (MeV 2 ) 


Figure 2: Kaon electromagnetic form factor in the space-like region. Solid 
line corresponds to the set of parameters fitted with m u = 350 MeV, dashed 
line corresponds to m u = 210MeV. The figure shows up the difference in 
the slope of the curves. 
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Figure 3: Pion electromagnetic form factor in the space-like region compared 
with experimental data |5Bj. Solid line corresponds to the set of parameters 
fitted with m u = 350 MeV, dashed line corresponds to m u = 210 MeV. 

why the difference in the kaon charge radius evaluated with these two sets of 
parameters is greater than it seems from fig. 1. Figure 3 presents the Pion 
electromagnetic form factor. 

In figure 4 we show the Kaon weak decay constant for the up constituent 
mass in the range 210MeK < m u < 350 MeV. We can see that f K decreases 
smoothly as m u increases. 

In order to investigate if a different choice of X q could provide a better 
fitting of the kaon weak decay constant in the m u = 350 MeV case, we also 
plot, in figure 5, the dependence of the kaon weak decay constant with \ q . 
We verify that it is still not possible to fit the Kaon weak decay constant 
with its experimental value for any value of X q . In fact, it is not possible to 
adjust the parameters in order to reproduce rn K = 494 MeV for \ q < 0.93. For 
0.93 < A q < 1.00 we have f K < 71.6 MeV, and as X q tends to infinity, the kaon 
weak decay constant tends to f K = 9A.6MeV , and the constituent strange 
mass, which also depends on X q , runs from rn s = QAAMeV to m s = 350 MeV. 
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120 



Figure 4: Kaon weak decay constant for different quark up constituent 
masses. 



k 

q 


Figure 5: Kaon weak decay constant as a function of X q for m u = 350 MeV. 
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6. Conclusions 


We investigated the role of ambiguities associated with finite surface terms 
of the type of differences between logarithmic divergencies in the LSM with 
fermions. We found in the pseudoscalar-axialvector sector an arbitrary c- 
number which cannot be absorbed by renormalization nor fixed through 
the Goldberger-Treiman relation or PCAC. In this process, however, nec¬ 
essary chiral and gauge invariant counterterms were introduced which link 
the pseudoscalar-axialvector sector with the scalar-vector one. The coun¬ 
terterms have a common coefficient which is fixed by the value assigned to 
the arbitrary c-number. If this number is chosen to be non-zero, it induces 
a symmetry breaking term in the scalar-vector sector. We must conclude 
therefore that the apriori arbitrary c-numbers are finally removed in order 
not to violate chiral symmetry in this sector. Further, we have obtained 
numerical results which show a good agreement between the model and ex¬ 
perimental data for the pion and kaon electromagnetic form factors, even 
for higher values of the quark up constituent mass. The kaon weak decay 
constant for rn u = 350 MeV, nevertheless, cannot be adjusted close to its 
experimental value in the range of the other free model parameters. 
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Appendix 


Consider first the following example. The surface-term integral 


vaf3 (jTl ) 


f d 4 p 

9 [iv a0 

n, PnPuPaP/3 1 

/ (2vr ) 4 

1 

to 

1 

to 

to 

-1 

CM 

1 

CM 

3 


(81) 
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does not depend on m 2 , i.e. 


dY, 


livafi 


dm 2 


= 2 ■ f d A p _ , 0 P^PuPaP /3 

% J (2tt) 4 [(p 2 *m 2 ) 3 (p 2 — m 2 ) 5 


= 0, 


(82) 


where we use the standard notation 

9fJ,ua/3 9^u9afS T 9^ia9v(5 “I - 9fj,f39va ■ 


Let us generalize Eq.(J82j). For that, note that in the space of d = 2a; 
dimensions one has 


n 1/42 M2n-lM2n _ On 

y^l/^ 2 ---At 2 ny • • • y — ^ 


r(o; + n) 

fM 


Therefore one obtains the identity 

tiuJPuiPi* ■ ■ ■ P» 2 „ _ RPT (A - u; - n) 


J d^p 1 - 


(p 2 — m 2 ) A 


(m 2 ) A ~ u ~ n 2 n r (A) 


9 /ll/J2.../12n 


(83) 


£Xj; n ) f j2oj 9llUJ.2-IJ.2n 

2 n T(A) J P (p 2 -m 2 ) A ~ 


Eq.( 


can be rewritten as 

9mH2---lJl2r, 


J d 2u p 


_ 2'T(v4) p Ml p M2 ■ ■ ■ p M2? 

(p 2 — m 2 ) A ~ n T(y4 — n) (p 2 — m 2 ) A 


= 0. 


(84) 


Integrating now with respect to m 2 we obtain the most general form of all 
possible surface terms 


qA- 2 ) = j [ 

‘->1/22—M2 n J 


d u p 

W) 1 


fl , MlM2---M2n 


2 n r(A-i) p /m p /J2 .. .p „ 2 


(p 2 — m 2 ) 71- ” -1 r(A — n — 1) (p 2 — m 2 ) 74-1 


(85) 
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